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REPORT  SUMMARY 


DIRECT  COOLING  OF  PROPULSION  DRIVES 
FOR  HIGH  POWER  DENSITY  AND  LOW  VOLUME 

Award  No.  N000 14  -01  -1-0634 


Background 

Thermal  management  is  presently  a  major  issue  in  the  power  electronics  field  due  to  the  exponential 
rise  in  power  densities,  and  hence,  increases  in  loss  density  and  device  temperatures.  The  ever  escalating 
thermal  management  problems  will  not  be  resolved  unless  power  semiconductor  devices  are  able  to 
operate  at  high  temperatures  (e.g.,  wide  band-gap  materials)  or  a  packaging/cooling  technology  is 
developed  that  extracts  effectively  enormous  amounts  of  waste  heat  from  the  electronics.  Recently,  spray 
cooling  has  gained  much  attention  in  various  sectors  because  of  its  ability  to  manage  anticipated  future 
heat  loads  (>  100  W/cm2). 

Objectives 

The  main  objective  of  this  grant  was  to  investigate  the  feasibility  of  using  spray  cooling  as  the 
thermal  management  technology  for  power  semiconductor  devices  used  in  power  converter  applications, 
in  particular,  naval  applications. 

Approach 

The  approach  followed  by  the  researchers  from  the  University  of  Arkansas  (UAF)  was  the  following: 

•  To  understand  the  basics  of  spray  cooling  modeling  involving  phase  change  since  previous 
research  efforts  have  been  mainly  experimental  with  little  research  done  on  a  comprehensive 
modeling  of  the  phenomenon. 

•  To  develop  a  power  packaging  methodology  for  power  semiconductor  devices  that  it  is 
compatible  with  spray  cooling.  The  research  team  felt  that:  (a)  spray  droplets  impacting  directly 
wire  bonds  and  bare  die  might  affect  adversely  the  reliability,  and  (b)  the  thermal  resistances 
between  the  bare  die  and  the  base  plate  must  be  reduced  to  maximize  the  benefits  of  spray 
cooling. 

•  To  develop  a  power  module  as  part  of  a  technology  demonstrator.  The  research  team  decided  to 
develop  an  H-bridge  topology  to  drive  a  dc  motor. 

•  To  perform  thermal  modeling  of  the  H-bridge  power  module  to  aid  in  the  design  stages. 

•  To  design  and  build  a  spray  box  for  the  developed  power  module  that  includes  not  only  a  module 
but  also  the  gate  drivers. 

•  To  test  the  technology  demonstrator  in  a  laboratory  environment. 

Insulated  Gate  Bipolar  Transistors  (IGBTs)  were  selected  as  the  power  semiconductor  device  of 
choice  because  of  their  availability  in  bare  die  form  from  manufacturers  at  reasonable  costs  and  the 
suitability  of  this  semiconductor  type  for  high  power  applications  up  to  2  MW  power  ratings. 

Conclusions 

The  conclusions  could  be  divided  into  two  groups;  namely,  those  related  to  spray  cooling  including 
theoretical  and  experimental  ones,  and  those  related  to  power  packaging  for  spray-cooled  power 
semiconductor  devices. 

Within  the  first  group,  the  researchers  established  that 

•  The  whole  process  of  modeling  phase  change  under  spray  cooling  is  a  difficult  task. 


•  Modeling  of  the  heat  transfer  mechanism  by  a  single  droplet  on  a  hot  surface  using  a  heat 
transfer  coefficient  and  a  CFD  package  has  limited  application  because  of  the  difficulty  in 
determining  the  heat  transfer  coefficient  under  spray  cooling  from  experimental  work. 

•  Multiphase  flow  modeling  is  a  viable  technique  to  understand  heat  removal  at  the  interfaces 
of  liquid  and  solid  by  a  vapor  bubble . 

•  The  “level  set”  method  is  one  of  the  best  methods  to  solve  multiphase  flow  in  spray-cooling 
systems. 

•  High  nozzle  flow  rates  are  he^>ful  when  all  other  conditions  are  the  same  because  of  the 
enhanced  conduction  mechanism  for  heat  removal.  However,  high-flow  nozzles  do  not 
always  produce  the  optimum  cooling  efficiency  since  they  may  hinder  the  phase-change 
process. 

•  Droplet  size,  drop  let  velocity  and  perhaps  other  factors  clearly  play  a  key  role  in  this  process. 
It  is  a  challenge  to  relate  theoretical  variables  such  as  droplet  size  and  droplet  velocity  to 
variables  that  are  used  in  the  design  of  spray-cooling  systems. 

•  Spray  cooling  is  much  better  than  liquid  cooling  of  a  standard  TO-220AB  package  containing 
IGBTs. 

Simply  spraying  a  conventional  power  module  does  not  result  in  the  optimum  cooling  efficiency. 
The  package  must  be  also  optimized.  Within  the  second  group  of  conclusions,  the  UAF  researchers 
accomplished  the  following: 

•  Eliminated  wired  bonds,  a  major  source  of  failure  in  power  modules,  by  replacing  them  with 
solder  balls. 

•  CTE  matched  the  package  materials  for  high  reliability. 

•  Reduced  the  thermal  resistance  between  the  heat  source  and  the  liquid  spray  significantly, 
and  thus  improving  overall  cooling  efficiency. 

•  Developed  a  power  packaging  technology  suitable  for  chip-scale  packages  so  power  density 
can  be  increased  further. 

•  The  initial  testing  of  the  spray-cooled  UAF  power  package  far  exceeded  the  performance  of 
the  liquid-  or  even  spray-cooled  conventional  package. 
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DIRECT  COOLING  OF  PROPULSION  DRIVES 


FOR  HIGH  POWER  DENSITY  AND  LOW  VOLUME 


I.  INTRODUCTION 

This  Final  Report  summarizes  the  research  work  and  conclusions  related  to  Award  N00014-01-1- 
0634  from  the  Office  of  Naval  Research  (ONR)  under  the  DoD  Experimental  Program  to  Stimulate 
Competitive  Research. 

The  research  team  from  the  University  of  Arkansas  at  Fayetteville  (UAF)  consisted  of  Drs.  J.  C. 
Baida,  F.  D.  Barlow,  R.  P.  Selvam  (Computational  Mechanics  Laboratory),  A.  Elshabini,  V.  Wang  and  K. 
C.  Burgers  from  the  Department  of  Electrical  Engineering  and  graduate  students  Mr.  Kiran  Vanam 
(doctoral  candidate),  Mr.  Jeremy  Junghans  (MS  candidate),  Ms.  Sandya  Bhaskar  (MS  candidate),  Mr. 
Suchit  Reddy  Damugatla  (MS  candidate)  and  Mrs.  Y.  Wei  (MS  candidate).  The  four  co-principal 
investigators  were  Dr.  Baida,  Dr.  Barlow,  Dr.  Selvam  and  Dr.  Elshabini.  Dr.  Wang  helped  with 
packaging  experimental  issues  and  Dr.  Burgers  designed  a  gate  driver  for  Insulated  Gate  Bipolar 
Transistors  (1GBT).  Ms.  Bhaskar  worked  on  theoretical  understanding  of  the  spray  cooling  phenomenon. 
Mr.  Junghans  developed  the  design  of  the  spray-cooling  system  and  Mr.  Vanam  the  packaging 
methodology  for  spray-cooled  power  modules.  Mr.  Damugatla  produced  FlowthermIM  thermal  models  of 
a  spray-cooled  H  bridge.  Mrs.  Wei  selected  the  cold  plate  to  compare  the  performance  of  a  packaged 
IGBT  when  cooled  by  the  cold  plate  versus  the  spray-cooling  system. 

Thermal  management  is  presently  a  major  issue  in  the  power  electronics  field  due  to  the  exponential 
rise  in  power  densities,  and  hence,  increases  in  loss  density  and  device  temperatures  (see  Azar  2002).  The 
ever  escalating  thermal  management  problems  will  not  be  resolved  unless  power  semiconductor  devices 
are  able  to  operate  at  high  temperatures  (e.g.,  wide  band-gap  materials)  or  a  packaging/cooling 
technology  that  extracts  effectively  enormous  amounts  of  waste  heat  from  the  electronics.  Recently, 
spray  cooling  has  gained  much  attention  in  various  sectors  because  of  its  ability  to  manage  anticipated 
fiiture  heat  loads  (>  100  W/cm2).  Thus,  the  main  objective  of  this  grant  was  to  investigate  the  feasibility 
of  using  spray  cooling  as  the  thermal  management  technology  for  power  semiconductor  devices  used  in 
power  converter  applications,  in  particular,  naval  applications.  From  the  initial  research  activities,  the 
research  team  established  the  following  enabling  objectives: 

•  Understand  the  basics  of  spray  cooling  modeling  involving  phase  change  since  previous  research 
efforts  have  been  mainly  experimental  with  little  research  done  on  a  comprehensive  modeling  of 
the  phenomenon . 

•  Develop  a  power  packaging  methodology  for  power  semiconductor  devices  that  it  is  compatible 
with  spray  cooling  (i.e.,  it  was  thought  that  spray  droplets  impacting  directly  wire  bonds  and  bare 
die  might  affect  adversely  the  reliability).  Previous  research  efforts  have  sprayed  the  base  plate 
of  the  power  module;  thus,  the  thermal  resistance  of  the  package  remained  unaltered  (see  Shaw 
(2002)).  The  research  team  felt  that  the  thermal  resistances  between  the  bare  die  and  the  base 
plate  must  be  reduced  to  maximize  the  benefits  of  spray  cooling.  IGBTs  were  selected  as  the 
power  semiconductor  device  of  choice  because  of  their  availability  in  bare  die  form  from 
manufacturers  at  reasonable  costs  and  the  suitability  of  this  semiconductor  type  for  high  power 
applications  up  to  2  MW  power  ratings. 

•  Develop  a  power  module  as  part  of  a  technology  demonstrator.  The  research  team  decided  to 
develop  an  H-bridge  topology  to  drive  a  dc  motor.  The  latter  was  chosen  because  of  its 
simplified  control  requirements  (as  opposed  to  those  for  an  induction  motor).  The  goal  was  to 


1 


evaluate  the  feasibility  of  spray  cooling  as  a  thermal  management  tool  for  power  semiconductor 
devices  independently  of  the  driven  electric  motor. 

•  Perform  thermal  modeling  of  the  developed  H-bridge  power  module  to  aid  in  the  design  stages. 

•  Design  and  build  a  spray  box  for  the  developed  power  module  that  includes  not  only  a  module 
but  also  the  gate  drivers  (which  were  designed  using  off-the-shelf  components). 

•  Test  the  technology  demonstrator  in  a  laboratory  environment. 

This  final  report  is  organized  as  follows:  Section  I  provides  a  general  overview  of  the  research  team 
and  the  project  objectives,  Section  II  presents  the  status  of  computer  modeling  of  the  spray  cooling 
phenomenon  during  the  period  of  performance  of  the  award.  Section  III  describes  the  philosophy 
developed  for  packaging  power  semiconductor  devices  subjected  to  spray  cooling.  Section  IV  briefly 
summarizes  thermal  modeling  of  the  H-bridge  module  using  Flowtherm™,  a  Computational  Fluid 
Dynamic  (CFD)  software  package  which  allows  simplified  simulation  of  the  spray  cooling  phenomenon, 
Section  V  describes  the  design  and  construction  of  a  spray  cooling  box  for  a  H-bridge  power  module 
intended  to  drive  a  dc  motor,  Section  VI  presents  the  experimental  results.  Section  VII  provides  the 
conclusions  of  the  research  work  enabled  by  this  ONR  research  grant,  and  finally.  Section  VIII  provides  a 
brief  summary  of  future  work  by  the  UAF  researchers. 
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II. 


COMPUTER  MODELING  OF  SPRAY  COOLING 


For  efficient  design  of  spray  cooling  systems,  understanding  of  the  spray  dynamics  and  droplet 
evaporative  cooling  phenomenon  via  numerical  simulation  is  very  important.  The  entire  physical 
processes  associated  with  this  phenomenon  are  very  complex.  The  processes  ranging  from  the  time  the 
fluid  comes  out  of  a  nozzle  orifice  as  spray  of  fine  droplets  that  will  impinge  on  a  hot  surface  and 
vaporize  to  the  time  that  the  vapor  is  removed  have  to  be  modeled  accurately.  There  is  not  much 
theoretical  work  done  to  theoretically  understand  the  spray  cooling  in  the  two-j)hase  regime  on  a  closed- 
loop  system  using  a  pure  single  component  medium  such  as  water  or  Fluorinert M  fluids  without  the  use  of 
air  as  atomizer  fluid.  Some  related  discussion  is  found  with  regard  to  fuel-air  spray  in  automobile  engine 
combustion  as  reported  by  Srignano  (1999).  DiMarzo  (1993)  modeled  the  evaporation  of  a  droplet  from  a 
hot  surface  for  different  applications.  Several  other  modeling  efforts  are  reported  in  the  literature  with 
regard  to  the  effect  of  the  impact  on  the  droplet,  associated  shape  change  mechanism,  and  other  issues. 
The  status  of  spray  cooling  was  surveyed  as  part  of  this  ONR  award  with  a  summary  given  below. 
Conclusions  and  future  work  on  theoretical  understanding  and  computer  modeling  of  spray  cooling  to  be 
addressed  in  a  followed  up  Phase  11  project  (ONR  Award  No.  N000 14-04- 1-0603)  are  given  in  Sections 
VII  and  VIII. 

11. 1  Effect  of  a  Single  Droplet  on  Heat  Transfer  Using  Heat  Transfer  Coefficient 

The  effect  of  a  single  droplet  on  heat  transfer  was  first  modeled  as  the  beginning  to  theoretically 
understand  the  spray  cooling  phenomenoa  The  interface  between  the  droplet  and  vapor  was  modeled 
using  a  moving  boundary  similar  to  the  work  reported  in  Gordon  (1999)  and  DiMarzo  (1993).  So  the 
amount  of  liquid  vaporized  at  the  interface  of  the  liquid-vapor  is  taken  away  from  the  droplet  at  each 
computing  time.  DiMarzo  measured  experimentally  the  convective  heat  coefficient  to  be  used  for 
modeling  purposes.  Their  modeling  assumed  that  the  vaporization  happened  at  the  interface  of  the  liquid- 
vapor.  Whereas  in  the  real  world,  depending  on  the  size  of  the  droplet,  evaporation  can  happen  through 
several  mechanisms  and  a  vapor  bubble  can  form  at  the  interface  with  the  liquid  or  at  an  interface 
involving  solid  and  liquid.  Other  models  were  investigated  since  the  reported  liquid-vapor  model  is  not 
representative  of  what  is  actual  happening  in  the  spray-cooling  phenomenoa 

11.2  Theoretical  Issues  in  Spray  Cooling 

There  are  several  theoretical  challenges  at  this  time  ranging  fom  the  efficient  design  of  jets  to 
convert  the  liquid  into  spray  and  finally  the  process  of  removing  the  heat  from  the  hot  surface.  There  is 
no  comprehensive  theoretical  model  available  at  this  time  for  spray  cooling  that  includes  phase  change. 
Hence,  there  are  no  CFD  packages  which  the  researchers  or  end  users  can  commercial  obtain  for  this 
application.  Complete  modeling  of  the  heat  removal  by  micro  jet  array  system  using  the  Finite  Difference 
Method  (FDM)  for  a  single  fluid  was  reported  by  Selvam  and  Elshabhini  (2000).  Availability  of  similar 
models  for  spray  cooling  should  be  a  valuable  tool  for  the  design  of  spray-cooling  systems.  Tryggvason 
et  al.  (2001)  gave  a  general  review  of  multiphase  flow  modeling  including  current  limitations  as  well  as 
specific  developments  in  the  front-tracking  method,  used  for  multiphase  flow  simulation.  They  discussed 
the  application  of  the  model  to  bubbly  flows,  atomization  of  liquid  for  combustion,  secondary  breakup  of 
drops,  primary  breakup  of  a  jet,  solidification  and  boiling  flows.  This  is  a  good  survey  paper  for  a 
beginner  in  this  area.  Dhir  (2001)  give  a  summary  of  their  numerical  simulation  for  pool  boiling  heat 
transfer  including  some  issues  related  to  microgravity  environments.  Particular  technical  issues  are 
specifically  addressed  in  the  next  few  sections. 

II.2.1  -  Flow  in  spray  nozzle:  Rahman  (2004)  analyzed  the  single -phase  flow  in  a  swirling  jet  using  a 
Finite  Element  Method  (FEM)  for  optimum  design.  The  numerical  simulation  predicts  the  jet  cone  angle 
and  pressure  drop  across  the  nozzle.  Further  work  is  needed  to  consider  the  formation  of  the  spray  and  its 
effect  in  the  design  of  jets. 
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11.2.2  -  Atomization  of  droplet:  Tryggvason  (2001)  g^ve  a  survey  on  computer  modeling  of  atomization 
of  droplets.  The  methods  were  classified  into  two  groups:  primary  breakup  and  secondary  breakup.  Few 
researches  including  Tryggvason’s  group  work  were  reported.  A  gpneral  survey  on  drop  and  spray 
formation  from  a  liquid  jet  was  provided  by  Lin  and  Reitz  (1998)  including  the  status  of  experimentation 
and  modeling. 

11.2.3  -  Droplet  and  surface  interaction:  Kizito  (2004)  used  the  Tryggvason’s  front  tracking  model  to 
study  the  splashing  droplet  on  a  hot  plate.  They  found  that  the  maximum  heat  is  transferred  on  the 
substrate  when  the  droplet  has  the  widest  extent  of  spread.  Pasandideh-Fard  (2001)  used  the  Volume  of 
Fraction  (VOF)  method  for  similar  work.  In  these  two  references  the  fluid  is  below  the  range  of  boiling 
and  they  both  used  a  droplet  having  a  2-mm  diameter. 

11.2.4  -  Liquid  and  vapor  interaction  or  bubble  cfynamics:  Extensive  computer  modeling  work  has 
been  done  in  this  area.  The  surveys  of  Tryggvason  (2001)  and  Sussman  (1994)  are  significant 
contributions.  Here,  surface  tension  and  gravity  have  an  impact  on  the  dynamics  of  the  bubble  in  the 
liquid.  The  shape  of  the  bubble  can  change  depending  upon  which  factor  dominates.  There  is  extensive 
knowledge  available  in  this  area.  This  is  part  of  the  boiling  flows  and  spray  cooling  when  vapor  grows  as 
a  bubble  and  moves  in  the  liquid. 

11.2.5  -  Liquid,  vapor  and  surface  interaction  on  a  hot  surface  or  boiling  flows :  The  phenomenon  of 
nucleation  during  boiling  and  the  growth  of  a  vapor  bubble  on  a  hot  surface  are  not  well  understood  and 
only  few  modeling  efforts  have  been  initiated  in  recent  years.  The  modeling  efforts  can  be  classified  into 
two  categories:  (1)  Multiphase  flow  modeling  considering  conduction,  convection  and  phase  change 
phenomena,  and  (2)  Simplified  procedures  where  only  conduction  and  phase  changes  are  considered. 
Welch  and  Wilson  (2000),  Son  and  Dhir  (1998)  and  Juric  and  Tryggvason  (1998)  have  computed  the  film 
boiling  stage  using  multiphase  flow  modeling  This  stage  was  modeled  earlier  even  though  it  occurs  after 
nucleate  boiling.  The  reason  might  be  that  the  density  ratio  is  far  smaller  than  that  one  in  nucleate  boiling 
situations.  Some  research  work  has  been  done  on  nucleate  boiling  and  will  be  reported  below. 

11. 2. 5.1  -  Film  Boiling :  Welch  and  Wilson  (2000)  used  the  VOF  method  to  model  film  boiling.  Juric  and 
Tryggvason  (1998)  used  the  front-tracking  method  to  study  film  boiling  in  2D  grids  while  Shin  and  Juric 
(2002)  used  it  for  3D  grids.  They  could  only  consider  cases  with  low  density  ratios.  The  numerical 
method  could  not  allow  them  to  apply  their  approach  for  high  density  ratios  that  occur  in  nucleate  boiling 
flows.  Son  and  Dhir  (1998)  used  the  level  set  method,  which  is  similar  to  the  method  in  Sussman  ( 1994), 
to  compute  film  boiling.  The  computed  results  were  in  reasonable  agreement  with  the  experimental 
results. 

11.2.5.2  -  Nucleate  Boiling :  Lee  and  Nydhal  (1989)  computed  the  single  bubble  growth  in  nucleation 
boiling  by  considering  only  the  liquid  region  and  assuming  the  bubble  to  be  spherical.  Then,  Son  (1999) 
applied  their  film  boiling  model  to  compute  bubble  growth  during  nucleation  boiling.  The  computed 
bubble  size  was  around  1  mm  and  compared  well  with  their  experimental  results.  The  bubble  growth 
during  nucleation  in  their  model  was  mainly  due  to  the  micro  region  model  which  was  based  on  a  ID 
model  close  to  the  dry  out  area  on  the  hot  wall.  For  the  macro  region  where  the  Navier-Stokes  equations 
were  used,  the  contribution  to  phase  change  was  very  small.  Bai  and  Fujita  (1999)  and  Yoon  (2001)  also 
computed  the  nucleated  pool  boiling  using  the  Lagrangian  method  for  tracking  the  interface  and 
considering  only  the  Iquid  region.  This  approach  was  similar  to  that  of  Lee  and  Nydhal  except  the 
meniscus  model  in  Bai  and  Fujita  was  close  to  Son  (1999).  Kunugi  (2001)  used  the  VOF  method  for 
boiling  and  condensation  problem  without  much  detail  on  the  numerical  implementation  or  the  problems 
that  were  considered. 
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II. 2. 5. 3-  Simplified  Models:  Di  Marzo  (1993)  modeled  a  single  water  droplet  on  a  hot  surface.  He  (2001) 
proposed  a  model  for  boiling  heat  transfer  based  on  a  numerical  macro-layer  model.  Here,  the  boiling 
curve  was  reproduced  numerically  by  determining  the  macro-layer  thickness.  They  explained  that  the 
evaporation  due  to  the  growth  of  vapor  stem  was  the  main  contribution  to  the  total  heat  flux.  They 
combined  this  model  with  a  conduction  model  for  transient  analysis.  Stephan  and  Hammer  (1994)  and 
Kern  and  Stephan  (2003)  used  micro  and  macro  regions  for  modeling  heat  transfer  during  nucleation 
boiling.  In  the  micro  region,  a  ID  complex  fourth-order  equation  was  solved  to  determine  the  thickness 
and  width  of  the  region  with  only  conduction  considered  in  the  macro  region.  They  compared  computed 
and  experimental  results. 


1 1.2.6  -  Spray  cooling:  From  the  above  discussions  one  can  conclude  that  several  different  research 
efforts  have  been  undertaken  to  understand  different  phases  of  boiling  and  most  of  them  were  for  pool 
boiling.  None  of  them  can  be  applied  directly  to  the  understanding  of  spray  cooling.  In  the  case  of  spray 
cooling  which  is  of  interest  to  electronic  cooling,  conduction  of  heat  from  the  hot  wall  into  the  liquid, 
nucleation  and  growth  of  a  vapor  bubble  on  the  hot  surface  due  to  phase  change,  convection  of  liquid  and 
vapor  bubble  in  the  liquid  film,  and  impact  of  a  spray  droplet  on  both  the  thin  liquid  film  and  a  vapor 
bubble ,  all  happen  simultaneously  as  shown  in  Figure  1 .  The  above  processes  happen  in  the  midst  of 
surface  tension  and  micro  or  macro  gravity  conditions.  Son  (1999)  could  grow  bubble  diameters  larger 

than  1.25  mm  and  moved  them  smoothly 
into  a  7.5-mm  pool  of  liquid  above.  Also, 
most  of  the  phase  change  happened  in  the 
micro  region. 

In  spray  cooling,  it  is  very  rare  to  have 
a  liquid  film  thickness  greater  than  a  couple 
of  millimeters.  Furthermore,  the  liquid 
film  on  the  hot  surface  will  be  in  the  range 
of  40  pm  to  200  pm  for  spray-cooling 
systems  designed  for  high  heat  fluxes. 
Usually  the  thickness  can  be  about  20%  of 
the  diameter  of  the  spray  as  observed  by 
Kizito  (2004)  in  their  experimental  and 
computer  modeling  work.  Lin  and 
Ponnappan  (2003)  reported  around  40  pm 
and  Rini  (2002)  reported  in  the  range  of 
200  pm. 

Theoretical  understanding  of  the 
nucleate  boiling  in  thin  films  in  the  range 
of  40  pm  to  200  pm  is  not  available  at  this 
time.  As  a  start  to  understand  the  spray 
cooling  phenomenon  in  a  thin  liquid  film, 
the  range  of  70  pm  will  be  considered  by 
the  UAF  researchers.  The  growth  of  vapor 
bubbles  and  their  impact  on  heat  removal 
on  the  hot  surface  will  be  investigated  in 
Phase  II  of  this  project  on  spray-cooling  systems  for  power  semiconductor  devices. 
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II.3  Methods  to  Solve  Multiphase  Flow 

The  major  classifications  of  methods  to  solve  multiphase  flow  are  the  Lagrangian  method,  Arbitrary- 
Lagrangian-Eulerian  (ALE)  method  and  Eulerian  method.  Some  of  the  applications  of  multiphase  flow 
modeling  using  these  methods  will  be  discussed  below  with  the  advantages  and  disadvantages  of  each 
method.  Tryggvason  (2001)  gave  a  comprehensive  survey  and  Scardovelli  and  Zaleski  (1999)  presented 
a  survey  on  VOF  with  discussions  on  numerical  techniques.  In  the  same  journal  volume  of  Tryggvason 
(2001),  several  other  review  papers  on  multiphase  flow  modeling  are  available  and  are  very  informative. 
The  equations  were  solved  using  either  FDM  or  FEM.  Of  these  two  methods,  the  FDM  is  the  most 
popular  one.  The  primary  numerical  procedure  used  by  the  UAF  researchers  is  either  FDM  or  FEM  on 
equal  grid  spacing. 

11.3.1  -  Lagrangian  method  This  method  is  fairly  complex  as  reported  by  Tryggvason  (2001).  The 
method  was  used  for  vapor  bubble  growth  modeling  by  Welch  (1998)  using  FEM  and  Lee  and  Nydhal 
(1989)  for  nucleate  boiling  using  FDM  on  a  body-fitted  grid  system. 

11.3.2  -  Arbitrary-Lagrangian-Eulerian  method  Yoon  (2001 )  applied  the  ALE  procedure  to  study  the 
bubble  growth  and  departure  from  the  heated  surface.  They  used  the  mesh  free  method  which  is  a  recent 
development  in  engineering  for  the  Eulerian  calculations.  The  implementation  may  be  complex  when 
compared  to  FDM  or  FEM. 

11.3.3  -  Eulerian  methods:  Front  capturing  methods,  the  separate  boundary-fitted  grid  method  for  each 
phase,  the  front  tracking  method  and  the  hybrid  of  front-tracking  and  front-capturing  methods  are 
available  at  this  time  for  multiphase  flow  calculations  as  reported  by  Tryggvason  (2001).  The  primary 
methods  used  in  front  capturing  approaches  are  the  Marker  and  Cell  (MAC)  method,  the  VOF  method, 
the  level  set  method  and  the  combination  of  VOF  and  level  set  methods.  The  separate  boundary- fitted 
grid  method  is  only  suitable  for  simple  geometries. 

Of  the  methods  mentioned  above,  the  VOF  method  and  level  set  methods  are  applied  extensively  in 
the  literature  for  multiphase  flow  modeling.  Discussions  on  VOF  and  level  set  methods  are  available  from 
Ferziger  and  Peric  (2002)  and  a  good  discussion  on  different  applications  of  the  level  set  method  is  given 
by  Osher  and  Fedkiw  (2003).  In  the  VOF  method,  sharp  differentials  have  to  be  calculated  between  the 
interfaces.  The  color  function  in  the  VOF  method  is  solved  using  a  time  dependent  convection  equation. 
This  is  similar  to  computing  shock  flow  in  compressible  flow  analysis. 

To  alleviate  the  problem  of  calculating  sharp  differentials  in  the  level  set  method,  a  distance  function 
is  introduced  which  will  be  varying  gradually  in  each  phase.  The  equation  to  be  solved  is  the  same  as  in 
the  VOF  method  except  the  use  of  a  distance  function  instead  of  the  color  function.  Once  the  distance 
function  is  solved  the  color  function  is  recovered  from  the  sign  of  the  distance  function.  The  level  set 
method  is  good  for  calculating  interfacial  curvature  accurately  but  mass  conservation  is  not  satisfied. 
Combining  with  the  VOF  method  this  error  is  eliminated  as  reported  by  Son  and  Hur  (2002).  Also,  it  is 
easy  to  handle  breaking  and  merging  of  an  interface  such  as  the  liquid-vapor  interface  in  the  level  set 
method.  In  addition,  the  implementation  of  3D  problems  is  a  straight  forward  extension.  Therefore,  the 
level  set  technique  was  selected  for  interface  capturing  in  the  work  forming  part  of  a  continuation  project 
(that  is,  ONR  Award  No.  N000 14-04- 1-0603). 
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III.  POWER  PACKAGING  COMPATIBLE  WITH  SPRAY  COOLING 


Reliability  may  be  compromised  when  graying  a  liquid  onto  fragile  wire  bonds,  which  are  a  major 
source  of  failure  for  packages  for  semiconductor  devices  (see  Ng  2000).  If  the  spray  can  be  kept  perfectly 
clean  this  may  be  possible,  however,  in  a  real  system  the  presence  of  impurities  in  even  very  small 
quantities  will  erode  the  surface  of  the  die  and  wire  bonds.  As  a  result,  the  die  must  be  protected.  In 
addition,  conventional  packages  were  not  developed  to  realize  the  full  potentials  of  spray  cooling.  These 
packages  contain  considerable  thermal  resistance  when  compared  to  the  low  thermal  resistance  to  the 
ambient  provided  by  spray  cooling. 

In  order  to  address  these  issues,  the  UAF  researchers  developed  a  novel  packaging  methodology  to 
replace  the  conventional  wire  bonds  and  allow  for  the  package  to  be  directly  spray  cooled  ( Vanam  2005). 
The  new  packaging  methodology  enables  an  important  reduction  of  the  thermal  resistance  between  the 
heat  source  and  the  sprayed  fluid,  thereby  greatly  improving  the  overall  cooling  efficiency.  The  UAF 
packaging  team  decided  first  to  build  a  test  vehicle  employing  a  thermal  test  die  to  demonstrate  both  the 
potential  of  spray  cooling  for  power  electronics  and  the  feasibility  of  the  proposed  packaging  technology. 
The  test  vehicle,  which  consists  of  a  thermal  test  die  with  built  in  temperature  sensors,  enabled  the 
researchers  to  easily  change  the  heat  load  and  understand  the  spray  cooling  phenomenon  under  different 
scenarios.  Finally,  the  packaging  methodology  was  implemented  in  an  H  bridge,  a  benchmark  circuit, 
using  IGBTs  from  International  Rectifiers  (IR).  The  thermal  capabilities  of  this  H-bridge  power  module 
were  demonstrated  by  driving  a  240- V,  1 0-hp  dc  motor.  A  patent  filing  is  under  preparation  for  this 
power  packaging  methodology  that  is  described  next. 


III.  1  ProposedPower  Packaging  Methodology 

The  proposed  base  package  shown  in  Figure  2  consists  of  a  heat  spreader  with  the  power 

semiconductor  device  (i.e.,  the  IGBT) 
Pottinir  romnound  flip  chipped  onto  a  substrate  such  as 
Heat  spreader  /  w  DBC  (Direct  Bond  Copper).  The  heat 

spreader  not  only  provides  an 
electrical  path  but  also  a  good  thermal 
path.  A  potting  compound  can  be 
used  for  encapsulation  of  the  power 
device.  The  top  and  bottom 

connections  of  the  power  device 
Figure  2.  Proposed  package  compatible  with  spray  cooling.  (IGBT)  with  the  heat  spreader  and  the 

DBC  are  created  using  a  high 

temperature  solder  and  a  low  temperature  solder,  respectively.  The  advantages  of  such  an  approach  are: 


•  Wire  bond  connections  are  eliminated 

•  Thermal  paths  are  improved 

•  High  heat  fluxes  in  excess  of  100  W/cm2  are  possible 

•  Double-side  cooling  is  possible 


Different  variations  of  the  base  package  are  possible  depending  upon  the  choice  of  the  heat  spreader 
material.  Some  choices  for  the  heat  spreader  material  are  diamond.  Metal  Matrix  Composites  (MMC), 
AIN,  and  CuMo.  In  the  case  of  diamond  or  AIN,  which  are  electrical  insulators,  modifications  are 
required  to  achieve  the  top  electrical  connection  (typically  the  collector  connection  for  an  IGBT).  The 
package  employing  diamond/AIN  as  the  heat  spreader  is  shown  in  Figure  3.  In  this  package,  the  kovar 
spacer  attached  between  the  heat  spreader  and  the  DBC  completes  the  top  electrical  connection. 
Appropriate  metals  are  deposited  to  achieve  solderable  connections  on  the  heat  spreader  and  kovar  spacer. 
This  approach,  apart  from  eliminating  the  wire  bond  connections,  also  reduces  the  electrical  resistances  of 
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the  package  interconnections.  This  reduction  allows  higher  current  capability  with  reduced  losses,  and 
hence  increases  power  density. 


Figure  4  shows  a  comparison  of  the  resistance 
bonds  and  solder  balls.  The  total  resistance  is  the 


Diamond  heat  spreader 


Potting  compound 


Kovar  leg 

E  -  Emitter 
G  —  Gate 
C  -  Collector 


Figure  3.  Diamond  as  a  heat  spreader  for  the 
proposed  package. 


for  the  emitter  contact  for  different  diameters  of  wire 
effective  resistance  seen  for  the  ‘n’  number  of  wire 
bonds/solder  balls  for  a  given  geometry.  Clearly 
the  resistance  offered  by  the  wire  bonds  is  far 
higher  than  that  of  the  solder  balls.  For  the  same 
diameter  of  1 0  mils,  the  resistance  of  wire  bonds  is 
120  times  greater  than  that  of  solder  balls.  Figure  5 
compares  the  total  resistance  of  the  wire-bonded 
module  and  the  flip-chip  package  (MMC  case). 
The  device  resistance  is  not  shown  as  it  is  assumed 
the  same  for  both  cases  and  is  excluded  in  the 
comparison.  It  shows  that  the  total  resistance  of  the 
wire-bonded  module  is  37  times  greater  than  the 
one  for  the  MMC  case. 


Diameter  (mil*) 


Wire  bonds 


E  // 

Solder  balk 


Die  area  =15.6  mm  x  15.6  mm 
Figure  4.  Comparison  of  resistance  between  wire  bonds  and  solder  balls. 


III.2  Fabrication  of  the  Test  Vehicle 

A  test  vehicle,  consisting  of  a  silicon- based  thermal  test  die  fabricated  by  the  researchers  and  a  CuMo 
heat  spreader  (35%  Cu  &  65%  Mo),  was  used  to  evaluate  the  proposed  packaging  methodology.  Table  1 
shows  the  thickness,  thermal  conductivity  K  and  Coefficient  of  Thermal  Expansion  (CTE)  for  different 
layers  in  the  package  that  is  closely  CTE  matched  to  maintain  the  mechanical  integrity.  The  thermal  test 
die  mimics  the  steady-state  losses  in  an  IGBT  and  has  embedded  sensors  to  determine  the  temperature 
across  the  die.  This  is  achieved  through  the  use  of  integrated  heaters  and  temperature  sensors  on  the 
silicon  die.  The  fabricated  test  die  measured  15  mm  x  15  mm  and  the  CuMo  heat  spreader  was  24  mm  x 
24  mm.  A  channel  was  milled  in  the  heat  spreader  to  accommodate  the  thermal  test  die  and  then  the  die 
was  attached  to  the  heat  spreader  using  a  high  temperature  solder.  Later,  the  die  was  flipped  and  attached 
onto  the  DBC  board  using  a  low  temperature  solder.  The  completed  test  vehicle  and  its  cross-sectional 
view  are  shown  in  Figures  6  and  7,  respectively. 
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Heat  spreader 


RMMC  3  2  ^\Rsold«-HT“  0.0004  pH 
^Sdder-Liv  0.25 

RSoldtr-LT  H7=  1  [if> 

Collector  Emitter 


Device  I  Leg 


DBC 

^total_MMC  =  3.25  J.l£2 
Wire  bonded  module 

RWue  Bonded- lOmil  =  1  - 1 


Figure  5.  Comparison  of  the  power  package  and  wire  bonded  module. 


Power  connector  for 
thermal  test  die 


CuMo  heat 


CuMo  heat  spreader 


Figure  6.  Completed  test  vehicle  (CuMo).  Figure  7.  Sectional  view  of  the  CuMo  package. 


Table  1  Thermal  and  mechanical  properties  of  various  layers. 


Interface 

Material 

Thickness 

(cm) 

K 

(W/cm-K) 

CTE 

(ppm/°C) 

Die 

Si 

0.045 

1.5 

4 

Solder 

95Pb  -  5Sn 

0.005 

0.23 

63Sn  -  T7RT- 

0.005 

25 

Heat 

Spreader 

65Mo—  35Cu 

0.3 

2.1 

7.2 

DBC 

AIN 

0.121 

1.7 

5 
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III.3  Fabrication  of  Flip-Chip  Power  Package  using  IGBTS  for  H-bridge  Implementation 

The  fabrication  of  the  power  package  can  be  divided  into  three  major  tasks: 

•  Processing  of  IGBT/diode  wafers  to  achieve  solderable  metallization 

•  Processing  of  the  substrate  (DBC) 

•  Assembly  of  the  package 

III.3.1  -  Processing  of  the  IGBT/diode  wafers 

Whole  wafers  of  the  IGBT  and  diode  die  were  obtained  from  International  Rectifier.  The 
1RG4PC30KD  IGBT  and  diode  die  were  selected  because  of  their  availability  in  wafer  form  and  were 
compatible  with  the  ratings  of  dc  motors  in  the  UAF  Energy  Conversion  Laboratory.  The  nominal  front 
metal  composition  (emitter  side)  for  IGBT  and  diode  (anode  side)  is  99%  A1  and  1%  Si,  which  is  suitable 
for  wire  bonding,  and  the  backside  metallization  is  Cr-Ni-Ag  that  can  be  soldered  to  a  substrate  using 
traditional  solders  such  as  Sn63.  Figures  8  and  9  show  the  diode  and  IGBT  wafers. 

As  aluminum  is  difficult  to  solder  to,  suitable  metallizations  are  required  for  soldering  the  front  metal 
during  the  flip  chip  process.  Therefore,  additional  thin  film  processing  of  the  obtained  wafers  was 
required.  The  major  additional  processing  steps  involved  are  listed  below: 

•  Dicing  of  the  wafer  for  processing  only  partial  portions  of  the  original  wafer 

•  Zincating 

•  Electroless  Nickel  (EN)  plating 

•  Dielectric  coating  (BCB/  Phtotoresist) 

•  Patterning  of  the  dielectric  layer  according  to  the  solder  mask  design 

An  EN  plating  process  was  chosen  to  plate  globally  the  exposed  metal  pads  as  it  (a)  minimizes  the 
number  of  process  steps  and  (b)  is  less  complicated.  However,  tight  control  of  the  plating  bath  is  required 
for  good  adhesion  of  the  nickel  onto  the  plating  surface.  This  means  that  a  nickel  plating  bath  along  with 
a  zincating  bath  were  setup  and  tested  on  several  test  dies  to  check  repeatability  and  process  optimization. 
The  bath  consists  of  25%  by  volume  of  the  zincate  concentrate  (NaOH,  zinc  oxide  and  sodium  gluconate) 
and  75%  of  distilled  water.  The  zincating  bath  chemically  removes  the  oxide  layer  on  top  of  the 
aluminum  and  at  the  same  time,  a  thin  layer  of  zinc  is  coated  on  the  aluminum  to  protect  it  from  forming  a 
native  oxide  layer.  As  the  zincated  sample  is  lowered  into  the  EN  plating  bath,  the  zinc  is  etched  away  by 
the  nickel  solution  and  plating  proceeds  onto  a  pure  clean  A1  surface.  The  EN  bath  consists  of  nickel 
sulfate  (part  A),  sodium  hypophosphite  and  ammonium  hydroxide  (part  B)  and  water  mixed  in  the  ratio  of 
1:2: 16.  The  bath  temperature  is  maintained  between  90-95  °C  by  a  ceramic  heater. 

Once  the  process  variables  were  optimized,  the  6inch  IGBT  wafer  and  4-inch  diode  wafer  were 
partially  diced  so  the  full  wafer  would  not  be  damaged  in  case  of  a  technical  problem.  Further  time 
delays  would  have  incurred  if  it  had  been  necessary  procuring  a  new  whole  wafer.  The  partial  wafers 
were  first  zincated  for  about  25  seconds  and  were  removed  after  the  grayish  appearance  of  zinc  on  the  A1 
exposed  pads.  Full  wafers  can  also  be  processed  in  this  manner,  however,  this  project  only  required  a 
small  number  of  dice  in  comparison  to  that  available  from  a  wafer,  and  the  wafer  cost  is  significant. 
Immediately  the  partial  wafers  were  placed  in  the  EN  bath  and  plated  for  about  10  minutes  to  yield  5  to  6 
microns  of  Ni.  The  EN  plated  wafers  are  then  thoroughly  rinsed  and  attached  to  a  5-inch  Si  carrier  wafer 
using  a  double  sided  thermal  release  tape.  This  method  allows  to  process  wafers  of  non-conventional 
sizes.  The  wafers  are  spin  coated  with  benzo-cyclo-butane,  which  is  a  dielectric  liquid,  and  patterned 
according  to  the  solder  mask  as  shown  in  Figures  10  and  1 1.  The  solder  mask  opening  size  was  evaluated 
from  6  to  9.5  mils,  prior  to  selection  of  the  9-mils  opening  since  this  size  resulted  in  the  best  reflow  of  the 
solder  balls  during  assembly.  The  patterned  partial  IGBT/diode  wafers  are  then  diced  to  yield  individual 
die.  Figure  12  shows  the  diced  IGBT  and  diode  pieces. 
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(a) 

Figure  8.  (a)  Whole  diode  wafer  (4-inch).  (b)  Bare  die  of  a  diode  and  (c)  Sectional  view  of  diode 


Figure  9.  (a)  Partially  diced  IGBT  wafer  (6  inch)  for  additional  processing,  (b)  Bare  die  of  an  IGBT 

and  (c)  Sectional  view  of  an  IGBT. 


III.3.2  -  Processing  of  the  substrate 

DBC  is  one  of  the  best  choices  for  high  current  carrying  requirements  (up  to  200  A).  This  is  due  to 
the  low  resistivity  of  copper  and  good  thermal  properties  of  DBC.  The  board  consisted  of  dielectric  (AIN 
or  AtO})  sandwiched  between  12-mils  of  pure  copper  metallization  on  either  side.  The  etching  mask  and 
the  schematic  for  the  H-bridge  board  are  shown  in  Fig.  13. 
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Figure  10.  AutoCAD  drawing  of  solder 
mask  for  diode. 


Figure  1 1 .  AutoCAD  drawing  of  solder  mask 
for  IGBT. 


A  liquid  photoimageable  solder  mask  (TA1YO  PSR-4000BN)  was  used  to  create  a  solder  mask  on  the 
DBC  boards.  The  solder  mask  opening  size  was  optimized  and  an  opening  of  9  mils  resulted  in  good 
adhesion  of  the  solder  balls  onto  the  substrate.  The  pitch  of  the  openings  was  selected  as  16  mils,  leaving 
a  gap  of  6  mils  between  two  adjacent  solder  balls  at  the  same  electrical  potential  Greater  isolation  values 
were  achieved  by  simply  spacing  the  solder  balls  an  appropriate  distance.  The  line  spacing  was  governed 
by  the  spacing  between  the  gate  and  emitter  pad  on  the  die  and  Figure  14  illustrates  the  line  spacing 
issues.  For  30-mils  gate  to  emitter  spacing  on  the  substrate  with  an  over  etch  of  4  mils,  a  row  of  solder 
ball  openings  were  etched  away  when  compared  to  the  15-mils  design  spacing  with  an  over  etch  of  2 
mils.  A  spacing  of  15  mils  was  selected  as  it  resulted  in  the  maximum  number  of  solder  ball  connections 
and  therefore  higher  current  handling  capability. 


EN  plated  openings  Dielectric  passivation 


Figure  12.  Diode  IGBT  with  passivation  and  solderable  metallization. 
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Figure  13.  Schematic  and  physical  layout  of  the  H-bndge  circuit. 


1 1 1.3.3  Assembly  of  the  power  package 

The  temperature  hierarchy  of  solder  is  the  key  for  integrating  the  whole  package  (i.e.  die,  heat 
spreader  and  substrate).  First,  a  set  of  four  diodes  and  four  1GBT  dies  were  reflowed  with  high 
temperature  95Pb-5Sn  solder  balls  (10  mils),  see  Figure  15.  Then  the  four  IGBT  and  diode  pairs  were  flip 
chipped  onto  the  DBC  substrate,  the  kovar  legs  (19-mils  thick)  were  also  placed  as  shown  in  the  Figure 
16.  The  kovar  legs  were  Ni  plated  (2  pm)  and  were  reflowed  along  with  the  die  using  95Pb-5Sn  solder 
performs.  The  final  stage  of  the  assembly  involves  reflowing  of  heat  spreader  (AIN)  onto  the  diode/IGBT 
pairs  using  a  low  temperature  solder  (Sn  63). 


6  mil  solder  mask  openings 


9  mil  solder  mask  openings 

Etched  row  of  solder  balls  on  the 
emitter  pad 


Figure  14.  Minimum  line  spacing  issues  on  the  substrate 


30  mil  spacing 


15  mil  spacing 


13 


Anode  pad 


Emitter  pad  Ga*Pad 


000^00 

ocoooo 


Figure  15.  Solder  balls  reflowed  on  the  Diode  and  IGBT  respectively 


The  side  of  AIN  onto  which  the  die  will  be  attached  was  sputtered  and  finally  plated  with  solderable 
metallization  (Ti/Cu/  Ni).  The  size  of  the  heat  spreader  was  1 1.58  mm  x  9.7  mm.  The  power  connectors 
(Cu  blocks)  were  also  soldered  at  this  stage.  Figure  17  shows  the  completed  H-bridge  power  module. 


Kovar  legs 

Figure  16.  Assembly  process  of  the  die  onto  the  substrate. 
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Heat  spreader  (AIN) 

Control  signal  inputs 

IiHrf/ 
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+  dc  bus  positive  terminal 

As\  U i  n.Witit  n  I.krmmo  1 

j^|4-  \1-  Motor  negative  terminal 

__ — __J  M  f  Motor  positive  terminal 

Figure  17.  Completed  10  hp  H-bridge  power  module. 
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rv.  THERMAL  ANALYSIS  OF  A  SPRAY-COOLED  POWER  MODULE  USING  FLOTHERM 


Thermal  modeling  of  the  developed  spray-cooled  H-bridge  power  module  was  performed  to  provide 
the  researchers  with  useful  information  regarding  the  potential  for  any  hot  spots.  This  effort  allowed  the 
researchers  to  optimize  the  design  before  building  prototypes.  Unfortunately,  there  was  no  commercial 
CFD  software  package  capable  of  simulating  the  two-phase  change  phenomenon  taking  place  under  spray 
cooling.  Although  commercial  packages  simulate  single -phase  cases,  they  can  be  used  to  estimate  the 
temperature  distributions  and  heat  paths  by  assuming  (a)  a  heat  transfer  coefficient  at  the  heated  surface 
impacted  by  the  droplets  (a  value  of  15,000  W/m2-K  was  used),  and  (b)  that  the  temperature  in  the 
chamber  (where  there  is  a  mixture  of  droplets  (fluid)  and  vapor  (gas))  is  constant  (this  temperature  was 
selected  as  40  °C  to  represent  the  worst  case  condition).  The  CFD  package  FlowthermrN1  was  used  to 
perform  the  thermal  simulations  of  the  developed  H-bridge  power  module  simply  because  it  was  available 
to  the  researchers. 

Using  the  device  datasheet  (i.e.,  part  IRG4BC30KD  from  IR),the  heat  losses  of  the  IGBT,  Pigbt,  can 
be  calculated  from  the  total  switching  energy  losses,  ESw ,  and  conduction  power  losses,  PCn ,  as  follows: 

PIGBT  =  ESwfs*  +  PCN  0) 

PcN  =  VcE(SAT)!  (2) 

where  fw  is  the  switching  frequency,  Vce(sat)  is  the  IGBT  collector-emitter  saturation  voltage,  Ic  is  the 
IGBT  collector  current,  and  D  is  the  duty  cycle. 

Similarly,  the  heat  losses  generated  by  the  freewheeling  diode,  P diode ,  are  calculated  from  the  diode 
reverse  recovery  losses,  Erec ,  and  conduction  losses,  Pcn_d>  as  follows: 

Pp/ODE  =  P  ree  f  nv  +  ?CN  _  D 

fcN.D  ~  —  D) 

where  VON  and  lR  are  the  on- state  voltage  and  on-state  current  flowing  through  the  freewheeling  diode, 
respectively. 

The  material  properties  of  the  different  blocks  in  the  power  module  are  the  following: 

Heat  Spreader 

■  Material:  AIN 

■  Dimensions:  1 1.58  mm  x  9.68  mm  x  0.625  mm 

■  Thermal  Conductivity:  170  W/m-K 

■  Specific  Heat:  725  J/kg-K 
Standoffs 

■  Material:  Kovar 

■  Dimensions:  2.0  mm  x  9.68  mm  x  0.7  mm 

■  Thermal  Conductivity:  16.3  W/m-K 

■  Specific  Heat:  439  J/kg-K 

Solder 

■  Material:  Sn63/Pb37 

■  Dimensions:  3.58  mm  x  4. 1 7  mm  x  0. 1  mm 

■  Thermal  Conductivity:  51  W/m-K 

■  Specific  Heat:  1500  J/kg-K 
IGBT  Bare  Die 

■  Material:  Silicon 

■  Dimensions:  3.58  mm  x  4.17mm  x  0.375  mm 

■  Thermal  Conductivity:  151  W/m-K 


(3) 

(4) 
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■  Specific  Heat:  750  J/kg-K 

Diode 

■  Material:  Silicon 

■  Dimensions:  3.43mm  x  2.54  mm  x  0.3  mm 

■  Thermal  Conductivity:  151  W/m-K 

■  Specific  Heat:  750  J/kg-K 
Solder  Balls 

■  Material:  Pb95/Sn5 

■  Radius  (each  one):  0. 125  mm 

■  Number  of  Solder  Balls:  8  x  36  =  288 

■  Thermal  Conductivity:  51  W/m-K 

■  Specific  Heat:  1500  J/kg-K 
DBC  Board 

Material:  Copper-AIN-Copper 
Copper 

■  Dimensions:33  mm  x  34.2  mm  x  0.3  mm 

■  Thermal  Conductivity:  385  W/m-K 

■  Specific  Heat:  385  J/kg-K 

■  Electrical  Resistivity:  3.2  x  10'x  Ohms-m 
AIN 

■  Dimensions:33  mm  x  34.2  mm  x  0.625  mm 

■  Thermal  Conductivity:  170  W/m-K 

■  Specific  Heat:  725  J/kg-K 

■  Electrical  Resistivity:  1.0x1 01 2  Ohms-m 


The  gaps  between  the  solder  balls  are  filled  with  an  underfill  having  a  thermal  conductivity  of  1 
W/m-K.  Figure  18  illustrates  a  side  view  of  the  H-bridge  power  module  implemented  in  Flowtherm  to 
determine  the  different  junction  temperatures,  7}.  By  using  the  above  formulae,  heat  losses  (inputs  on  the 
Flowtherm  simulations)  were  calculated  for  different  switching  frequencies  and  collector  currents. 
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Figure  19  shows  a  typical  Flowtherm  output  for  the  case  of  a  collector  current  of  18  A  and  a 
switching  frequency  of  6  kHz.  The  motor  is  rotating  in  only  one  direction  so  two  IGBTs  were  operating 
and  their  corresponding  freewheeling  diodes  operated  when  these  IGBTs  were  turned  off.  As  expected 
the  heat  losses  of  the  IGBTs  were  much  greater  than  those  of  the  freewheeling  diodes. 

Table  2  summarizes  the  results  for  switching  frequencies  of  6,  9,  12  and  15  kHz  as  well  as  collector 
currents  of  10,  16,  18,  20  and  30  A.  It  is  noted  that  16  A  is  the  maximum  continuous  collector  current  of 
the  selected  IRG4BC30KD  1GBT  for  a  case  temperature  of  100  °C.  Table  2  highlights  those  operating 
conditions  where  a  junction  temperature  of  125  °C  was  not  exceeded.  The  maximum  junction 
temperature  specified  by  the  manufacturer  is  150°C  but  the  research  team  decided  to  adopt  the  limit  of 
125  °C  due  to  safety  reasons.  For  the  IR-packaged  IGBT,  the  maximum  switching  frequency  is 
approximately  9.2  kHz  for  a  collector  current  of  16  A  a  junction  temperature  of  125  °C,  an  ambient 
temperature  of  40°C  and  other  data  (e.g.,  thermal  resistances)  from  its  datasheet 


Figure  19.  H-bridge  temperature  distribution  for  a  collector  current  of  18  A  at  6  kHz. 


Table  2.  Summary  of  the  Flowtherm-based  thermal  simulations. 
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The  continuous  collector  current  specified  by  the  manufacturer  depends  among  other  factors  on  the 
thermal  resistances  of  the  package,  a  TO-220AB  package  in  this  case.  Spray  cooling,  by  being  more 
effective  in  extracting  heat  losses,  may  offer  the  possibility  of  increasing  either  the  switching  frequency  or 
the  current  rating  of  the  semiconductor  device.  From  the  table,  spray  cooling  could  enable  operation  at 
the  higher  switching  frequency  of  15  kHz  for  the  rated  current  of  16  A  and  a  junction  temperature  of  125 
°C;  a  significant  increase  in  switching  frequency.  Alternatively,  spray  cooling  could  offer  the  possibility 
of  operating  the  IGBT  at  a  collector  current  of  18  A  for  a  switching  frequency  of  9  kHz  -  the  frequency 
limit  for  the  TO-220AB  package. 


19 


V.  DESIGN  OF  A  SPRAY  COOLING  BOX  FOR  A  FOUR-SWITCH  POWER  MODULE 

An  innovative  spray-cooling  system  illustrated  in  Fig.  20  was  developed  by  the  researchers  in  order 
to  demonstrate  the  heat-removal  capabilities  of  spray  cooling.  The  system  was  designed  to  cool  the  H- 
bridge  power  module  whose  design  and  packaging  methodology  were  addressed  in  Section  III  above. 
The  H  bridge  is  intended  for  driving  a  240- V,  10-hp  dc  motor  load. 

In  this  innovative  design,  the  spray  box 
contains  most  of  the  system  components 
allowing  for  a  compact  realization  in  order  to 
maximize  power  and  mass  budgets.  The  spray 
box  was  made  out  of  aluminum  with  an 
acrylic  airtight  lid.  The  fluid  is  contained  at 
the  bottom  of  the  spray  box  and  pumped  to  the 
nozzle  manifold  using  a  pump.  A  closed  loop 
was  formed  by  the  fluid  reserve  at  the  bottom 
of  the  box,  the  pump,  the  nozzle  manifold  and 
the  box  walls  acting  as  a  condenser. 

The  outer  dimensions  of  the  spray  box 
were  6.5”  x  6”  x  6”  (HxWxD)  with  the  H 
bridge  and  its  control  board  mounted  to  the  lid 
to  allow  accessibility.  The  motor  control 
circuitry  was  fabricated  on  a  FR-4  board  and 
placed  below  the  DBC  board  inside  the  spray 
box  as  illustrated  in  Figure  20.  The  motor 
control  board  has  a  1.2”  x  1.2”  window  (see 
Figure  21)  to  allow  for  the  spray  from  the 
nozzle  manifold  to  reach  the  H  bridge  located  above  and  facing  downwards.  The  manifold  shown  in 
Figure  22  was  designed  to  allow  one  nozzle  for  each  heat  spreader  (i.e.,  four  nozzles  in  total)  placed  at  the 
optimal  distance  below  the  H  bridge  (see  Wei  2004). 


Figure  2 1 .  Gate  driver  board  for  the  H  bridge. 
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H  bridge 


Figure  20.  Graphical  representation  of  the  spray  box. 


The  nozzles  were  type  1/8HH-1  from 
Spraying  Systems  Co.  (www.spray.com)  having  a 
nominal  orifice  diameter  of  0.031”,  a  length  of 
7/8”,  minimum  pressure  of  10  PSI,  and  maximum 
pressure  of  150  PSI.  A  single  nozzle  is  observed 
on  the  top  right  side  of  Figure  22.  The  airtight  lid 
also  contains  the  required  power  and  fiber  optic 
inputs  for  the  electronics.  IGBT  gate  signals  were 
optically  isolated  from  the  dc  motor  controller.  An 
observation  window  was  located  on  one  side  of  the 
spray  box.  The  window,  not  necessary  in  a 
commercial  product,  added  some  height  to  the 
spray  box.  Figure  23  shows  the  spray  box  housing 
with  it’s  observation  window. 

In  a  closed-loop  system,  the  vapor  that  results 
when  the  sprayed  fluid  in  the  form  of  droplets 
impacts  the  heated  surface  (or  device)  and 
undergoes  a  phase  change  must  be  recovered  by  a 
condenser.  Traditionally  the  two  tasks  of 
condensing  the  liquid  and  removing  the  heat  are  accomplished  by  adding  a  separate  condensation  coil  and 
heat  exchanger  to  the  spray  system.  This  methodology  has  proven  to  be  effective  but  significantly 
increases  the  overall  system  size,  thus  reducing  the  overall  system  power  density. 

The  proposed  spray  box  achieves  both  tasks  without  the  penalty  of  increased  system  size.  This  was 
accomplished  by  avoiding  the  use  of  a  dedicated  heat  exchanger  and  condensation  coil,  using  instead  two 
%”  cooling  channels  that  were  imbedded  into  the  Vi”  walls  of  the  spray  box  (not  shown  in  Figure  20). 

The  channels  allowed  chilled  water  to  be 
connected  to  the  spray  box  and  circulated 
throughout  the  aluminum  walls  in  a  fashion  very 
similar  to  that  of  a  cold  plate.  The  cooled  walls 
were  then  used  to  subcool  the  fluid  contained  at 
the  bottom  of  the  spray  box  while  condensing 
the  liquid  vapors  present  inside  the  box.  The 
unique  design  of  the  spray  box,  in  particular  the 
walls,  allows  the  system  to  maximize  heat 
removal  while  also  minimizing  the  size  and 
weight  of  the  spray-cooling  system. 

This  system  realization  is  different  from 
other  researchers  where  the  spray-cooling 
system  and  the  system  containing  the  heat 
sources  (the  H  bridge)  are  “packaged” 
separately  yielding  solutions  requiring  higher 
mass  and  volume  budgets  in  addition  to  the 
additional  thermal  resistance  of  the  separate 
package.  The  UAF  approach  consists  of 
designing  a  spray  box  that  becomes  the 
“package”  for  the  system  containing  the  heat 
sources. 
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Figure  23.  Completed  spray  box  housing. 


Before  designing  the  H-bridge  spray  system,  the  researchers  designed  and  developed  a  test  bench  to 
evaluate  the  feasibility  of  spray  cooling  for  power  modules.  The  Second  Annual  Report  covered  in  detail 
this  test  bench  (see  also  Wei  2004,  Ortiz  1999  and  Liu  2000).  The  main  purposes  of  the  test  bench  were 
to  allow  computation  of  key  parameters  such  as  heat  flux  (q”)  and  device  junction  temperature,  and  to 
optimize  the  spray  cooling  system  before  designing  the  final  system  consisting  of  a  spray-cooled  power 
module  with  real  power  semiconductor  devic  es. 

The  two  main  elements  of  the  spray-cooling  test  bench  are  the  test  bed  and  the  test  vehicle.  The  test 
bed  is  based  on  a  system  originally  used  to  spray  VME  cards.  The  system  was  modified  to  handle  higher 
power,  the  use  of  conductive  fluids  and  to  allow  quick  changes  in  the  experimental  setup.  A  single  spray 
nozzle  is  attached  to  a  linear  motor  allowing  simple  adjustments  of  the  spray  distance.  A  variable  speed 
motor  drives  the  pump  for  the  spray  loop.  This  allows  the  flow  rate  through  the  spray  nozzle  to  quickly 
be  changed  by  altering  the  speed  of  the  driving  motor.  A  heater/chiller  and  heat  exchanger  were  added  to 
the  spray  loop  enabling  precise  control  of  the  fluid  temperature.  Previous  work  by  the  UAF  researchers 
(see  Wei  2004)  indicated  that  the  ability  to  control  the  temperature  of  the  sprayed  fluid  is  significant  to 
achieve  optimal  heat  removal.  The  test  vehicle  consists  of  thermal  test  die  having  the  same  area  as  the 
selected  IGBT  and  mimicking  its  heat  losses. 
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VI.  EXPERIMENTAL  RESULTS 


The  experimental  results  have  been  divided  into  two  main  groups.  The  first  group  summarizes  the 
main  spray  cooling  results  that  were  obtained  using  (a)  a  thermal  test  die  and  the  test  bench,  and  (T^the 
spray  box  described  in  Section  V.  The  goal  of  the  test  bench  experimentation  was  to  evaluate  different 
type  of  nozzles,  fluids,  pressures  and  flow  rates.  The  advantage  of  spray  cooling  is  illustrated  by 
comparing  experimental  results  when  an  IGBT  packaged  by  the  manufacturer  is  cooled  by  fi)  the  spray¬ 
cooling  system  of  Section  V  and  (ii)  a  cold  plate.  The  second  group  corresponds  to  experimental  results 
related  to  spray  cooling  of  IGBT  bare  die  packaged  as  per  the  methodology  given  in  Section  III.  The 
main  goal  here  is  to  emphasize  that  the  final  thermal  management  solution  must  also  take  into  account  the 
packaging  aspects  if  one  is  to  fully  realize  the  potentials  of  spray  cooling. 


VI.l  Experimental  Results:  Test  Bench 

Previous  results  that  evaluated  several  spray  nozzles  from  Spraying  Systems  enabled  the  researchers 
to  select  nozzle  TG-1  since  it  yielded  the  best  results  for  this  application  when  using  a  mixture  of 
Fluorinert™  FC72  and  FC84  (see  the  Second  Annual  Report  and  Wei  2004).  The  saturation  temperature 
of  FC72  and  FC84  are  56  °C  and  80  °C,  respectively.  The  saturation  temperature  of  the  mixture  is  not 
known  but  the  boiling  point  of  the  mixture  tends  to  be  nearer  to  the  fluid  which  is  more  volatile  or  has 
higher  vapor  pressure  according  t>  the  manufacturer  3M;  that  is,  FC72  (232  \s.  79  torr).  The  specific 
heat  of  both  fluids  is  1.05  (J/g-°C),  the  latent  heat  and  density  of  FC-72  are  21  (cal/g)  and  1.68  (g/cnr  ), 
the  latent  heat  and  density  of  FC-84  are  19  (cal/g)  and  1.73  (g/cm3).  The  distance  between  the  spray 
nozzle  and  the  thermal  test  heater  was  varied  along  with  spray  pressure,  flow  rate  and  fluid  temperature. 
The  heat  flux  (q”)  was  calculated  from  the  observed  data  (see  Mudawar  2001)  as  follows: 


where  Pd  is  the  power  dissipated  by  the  thermal  test  heater  and  Ath  is  the  cross  area  of  the  heater  (i.e.,  15 
mm  x  15  mm).  In  this  case  the  sprayed  area  is  equal  to  that  of  the  test  heater.  Figure  24  illustrates  a  test 
heater  fabricated  by  the  UAF  researchers.  The  power  supplied  to  the  test  heater  was  slowly  increased 
until  it  reached  the  maximum  temperature  of  150  °C.  A  heater/chiller  and  heat  exchanger  were  used  to 
control  the  temperature  of  the  dielectric  fluid  being  sprayed  and  thermocouples  were  used  to  measure  the 
fluid  temperature.  Experiments  verified  that  heat  flux  is  maximized  for  given  nozzle  and  test  heater  when 
the  square  thermal  test  heater  was  circumscribed  by  the  circular  spray  pattern.  Thus,  the  optimal  distance 
is  a  function  of  the  spray  angle  and  heater  dimensions  (see  Wei  2004);  in  particular: 


0.5  wth 

sp  tan(O.50ip) 


(6) 


where  dv  is  the  optimal  spray  distance,  w,h  the  width  of  the  square  thermal  test  heater  and  0sp  is  the  nozzle 
spray  angle . 

Figure  25  illustrates  the  effect  of  the  fluid  temperature  upon  the  heat  flux.  Lower  fluid  temperatures 
(increased  subcooling)  resulted  in  higher  heat  fluxes.  This  may  indicate  that  in  addition  to  vaporization  of 
the  Fluorinert1™,  heat  removal  through  conduction  is  an  important  component  of  the  spray  cooling 
observed  in  the  test  bench.  The  maximum  heat  flux  q”  was  near  140  W/cm“  at  atmospheric  pressure  and 
a  flow  rate  of  0.085  gpm  for  the  TG 1  nozzle  that  has  a  0.036”  (36  mils)  orifice. 

Experiments  were  also  conducted  to  determine  the  role  of  spray  cooling  in  the  overall  heat  removal  of 
the  system.  The  total  heat  flux  of  the  system  (q”)  can  be  expressed  as  the  following: 


q”  =  q”s+q”c 


(7) 
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Figure  24.  Individual  test  heater  fabricated  at  UAF. 


where  qfs  is  the  heat  flux  associated  with 
spray  cooling  while  q”c  is  the  heat  flux 
due  to  conduction  losses  through  the 
backside  of  the  silicon  test  heater,  solder 
and  substrate.  Two  packages  were 
designed  and  built  h  order  to  determine 
the  significance  of  the  conduction  losses. 
The  first  package  utilized  materials 
chosen  to  maximize  the  conduction 
losses.  Highly  conductive  AIN  DBC  was 
chosen  as  the  substrate  and  95Pb-5Sn 
solder  as  the  die  attach  material.  The 
second  package  was  designed  to 
minimize  the  conduction  losses.  FR-4 
was  chosen  as  the  substrate  material  and 
non-thermally  conductive  adhesive  was 
used  to  attach  the  test  heater.  Both 
packages  were  spray  cooled  and  their 


total  heat  fluxes  are  compared  in  Figure  26.  The  results  indicated  that  just  as  high  total  heat  fluxes  were 
achieved  with  the  package  possessing  minimal  backside  conduction  as  with  the  highly  conductive 
package.  These  results  support  the  assumption  that  the  total  system  heat  flux  is  predominantly  a  function 
of  the  heat  removed  as  a  result  of  spray  cooling  (q”  ~  q”s). 


Similar  experiments  were 
conducted  leaving  fluid 
temperature  constant  (tf=  20  °C) 
and  varying  the  flow  rate.  The 
results  shown  in  Figure  27 
indicate  that  increasing  the  flow 
rate  led  to  higher  heat  fluxes  while 
the  other  parameters  remained 
constant.  These  experiments  were 
performed  using  a  single  type  of 
spray  nozzle. 

Similar  experiments  were 
performed  where  two  similar 
nozzles  having  slightly  different 
parameters  were  compared.  The 
experiments  compared  a  TG-1 
nozzle  (36  mil  orifice  diameter) 
and  a  TG-2  nozzle  (47  mil  orifice 
diameter).  The  results  given  in 
Figure  28  indicate  that  although 
the  TG-2  has  a  larger  capacity, 
thus  allowing  higher  flow  rates, 
nozzle  TG-1  results  in  better  heat 

fluxes.  This  indicates  that  the  specific  nozzle  parameters  such  as  orifice  diameter,  droplet  size  and  droplet 
velocity  play  a  significant  role  in  determining  the  amount  of  heat  flux  achieved  by  the  spray  system. 
These  results  also  indicate  that  conduction  is  an  important  component  of  the  arrangement  in  the  test 
bench.  As  mentioned  in  Section  II  above,  there  is  no  comprehensive  theoretical  model  available  at  this 
time  for  the  spray  cooling  phenomenon  that  can  allow  relating  parameters  like  flow  rate,  nozzle  orifice 
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Junction  Temperature  (°C) 

Figure  25.  fcttect  of  fluid  temperature  on  heat  fluxes  as 
function  of  heater  temperature. 
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q"  vs  Junction  Temperature 


— 0.19  gpm,  Tf=10,  AIN  0.26  gpm,  Tf=10,  AIN 

-A—  0.19  gpm,  Tf=10,  FR4  0.26  gpm,  Tf=10,  FR4 


Figure.  26  Heat  flux  comparison  of  highly  conductive 
(AIN)  package  and  poorly  conductive  (FR-4)  package. 


Figure  27.  Effect  of  flow  rate  on  heat  fluxes  as 
function  of  heater  temperature. 


and  others  into  parameters  like  droplet 
size,  droplet  velocity  and  thin  film 
thickness  that  are  used  for  computer 
modeling.  Making  such  a  connection 
will  be  a  goal  of  the  next  phase  of  this 
research  project. 

Power  Packaging:  Spray-Cooled  Test 

Vehicle 

The  feasibility  of  the  whole 
packaging  methodology  was  first 
investigated  using  a  low-cost  test 
vehicle  before  apply  ing  the  proposed 
packaging  techniques  using  1GBT 
bare  die  (see  Section  III.2  above). 
The  test  vehicle  consisted  of  a  thermal 
test  die  and  a  heat  spreader.  The 
thermal  test  die  mimicked  the  losses 
in  the  IGBT  and  had  embedded  sensors  for 
determining  the  temperature  across  the  die. 
Two  packaging  prototypes  were  used  for 
the  testing,  one  constructed  with  a  diamond 
heat  spreader  and  the  other  with  CuMo 
heat  spreader  (35%  Cu  and  65%  Mo). 
Each  package  was  sprayed  with  a  mixture 
of  50%  FC-72  and  50%  FC-84  diefectric 
fluids  at  standard  atmospheric  pressure,  a 
flow  rate  of  0.19  gpm  and  a  fluid 
temperature  of  10  °C.  The  experimental 
results  are  displayed  in  Figure  29.  Both 
packages  resulted  in  maximum  heat  fluxes 
well  over  100  W/cm:  for  a  maximum 
temperature  of  120  °C,  with  the  diamond 
package  exceeding  200  W/cm2.  These  heat 
fluxes  were  calculated  using  die  area  as 
opposed  to  sprayed  area  (the  latter  is  2.56 
larger).  The  setup  successfully 

demonstrated  the  feasibility  of  the 
proposed  spray  cooling  and  packaging 
methodology.  It  was  inferred  from  the  test 
bench  results  that  ai  atomized  fluid  that 


forms  very  fine  droplets  of  the  liquid 
produces  superior  results  when  compared  to  that  of  a  high  pressure  “fire  hose”  type  of  spray.  As  a  result 
the  UAF  researchers  selected  for  the  spray  box  (see  Figure  23)  the  nozzle  type  1/8  HH-1  that  produced  a 
finer  spray  than  the  nozzles  used  in  the  test  bench  experiments. 
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VI.2  Experimental  Results:  IR-Packaged  IGBT 

Bearing  in  mind  that  the  main  objective  of  this  grant  was  to  investigate  the  feasibility  of  using  spray 
cooling  as  the  thermal  management  technology  for  power  semiconductor  devices  used  in  power  converter 
applications,  in  particular,  naval  applications,  the  researchers  setup  a  second  test  bench  to  drives  a  240- V, 
10-hp  dc  motor.  Figure  30  is  a  photograph  of  the  test  bench.  A  dc  power  supply  (rated  0-300  V  and  50 
A)  is  not  shown  in  the  figure. 


q"  vs  Junction  Temperature 


— *—  TG-2,  0.38gpm,  Tf  =  10  TG-2,  0.59gpm.  Tf  =  10 

— * —  TG-1,  0.19gpm,  Tf  =  10-S-TG-1.  0.26gpm,  Tf  =  10 

Figure  28.  Effect  of  nozzle  parameters  on  heat  flux. 


Laboratory  experiments 
were  conducted  on  a 
commercial  TO-220AB 
package  of  the  IGBT 
IRG4PC30KD  using: 

•  The  spray  box  (see 
Figure  23),  and 

•  A  standard  cold  plate 
(see  Figure  31).  The 
IGBT  was  attached  to  the 
cold  plate  via  a  thermal 
pad  to  provide  electrical 
isolation  and  improve  the 
thermal  path  (due  to 
roughness  of  the  cold 
plate  surface). 


q"  vs  Junction  Temperature 


-  Diamond.  50/50  mix,  0.19  gpm,  Tf=10 


•  CuMo,  50/50  mix.  0.19  gpm,  Tf=10 


Figure  29.  Heat  fluxes  when  using  diamond  and  CuMo  heat 
spreaders  (die  area). 


For  the  sake  of 
simplicity,  only  one  switch  was  used  to 
run  the  dc  motor  and  the  temperature  of 
the  chiller  water  was  maintained  at 
around  20  °C  for  both  the  spray-cooling 
and  cold-plate  experiments  (the  chiller 
is  shown  on  the  right  side  of  Figure  30). 
A  mixture  of  50%  FC-72  and  50%  FC- 
84  dielectric  fluids  was  used  for  the 
spray-cooling  experiments. 

Tables  3  and  4  depict  the 
experimental  results.  The  collector 
current,  case  temperature  and  spray 
temperature  were  measured.  Power 
losses  were  calculated  as  mentioned  in 
Section  IV  using  switching  losses  and 
lowest  on-state  voltage  from  the  IGBT 
datasheet.  For  the  spray-cooling 
experiments  the  case  temperature  refers 
to  that  temperature  measured  on  the 
other  side  of  the  surface  where  the 
liquid  droplets  impinge.  That  is,  the 
side  of  the  base  plate  where  the  bare  die 
is  attached  (if  one  looks  at  the  package 
with  the  base  plate  facing  down  one 
could  see  two  areas  having  “half  circle” 
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shape  on  the  side  of  the  package  where  the  base  plate  is  visible).  The  junction  temperature  was  calculated 
using  the  measured  case  temperature  and  the  case-junction  thermal  resistance  from  the  IGBT  datasheet. 
All  calculated  values  were  rounded  off  to  their  nearest  integer  values  and  tie  operating  conditions 
common  to  both  cases  are  highlighted  in  yellow.  The  cold  plate  measurements  were  stopped  at  15  A 
because  the  measured  case  temperature  was  above  100  °C  indicating  that  the  junction  temperature  was 
higher  than  its  maximum  limit  of  1 50  °C  (as  it  was  later  calculated  after  completing  the  measurements). 
The  spray-cooling  measurements  were  stopped  at  19.5  A  because  of  a  limitation  of  the  dynamometer 
available  in  the  Energy  Conversion  Laboratory. 

Table  3.  Spray-cooling  data  for  IRG4BC30KD  IGBT  (1R  packaged)  at  9.2  kHz. 


lc(A) 

T(case)  (°C) 

T(junction) 

<°C) 

T(fluid  spray) 
(°C) 

Pigbt(W) 

0 

21.0 

21.0 

19.2 

0.0 

1 

23.6 

39 

19.3 

13 

10.4 

26.8 

59 

19.9 

15.2 

31.7 

79 

20.5 

40 

16.5 

33.2 

85 

20.9 

43 

19.5 

36.3 

97 

21.5 

51 

Table  4.  Cold- plate  data  for  the  IRG4BC30KD  IGBT  (IR  packaged)  at  9.2  kHz 


lc(A) 

T(case)  (°C) 

Tflunction) 

(°C) 

Pigbt(W) 

~ 0~ 

22.0 

22 

~s~ 

4.5 

48.1 

62 

12 

10 

81.9 

113 

26 

15 

118.1 

165 

39 

(a) 

Figure  30.  (a)  Experimental  setup  for  spray  cooling  of  power  devices,  (b)  Cooling  chamber. 
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IGBT  (underneath  the  board) 


Figure  31  Experunental  setup  for  liquid  cooled  cold  plate. 

Nevertheless,  a  comparison  of  the  measured  case  temperature  clearly  shows  the  advantages  of  spray 
cooling;  it  was  only  36.3  °C  at  a  collector  current  of  19.5  A  and  a  switching  frequency  of  9.2  kHz.  The 
above  operating  condition  would  definitely  be  outside  the  manufacturer  envelope  when  cooled  by  a  heat 
sink  or  a  cold  plate.  The  measured  case  temperature  for  the  cold -plate  case  is  around  272%  higher  than 
the  one  for  the  spray-cooling  case  for  approximately  15  A.  This  is  a  very  important  difference  that  could 
be  used  to  either  increase  the  switching  frequency  or  the  collector  current  of  the  IGBT.  The  reader  is 
reminded  that  the  maximum  continuous  collector  currents  specified  by  the  manufacturer  for  the 
IRG4PC30KD  are  16  A  and  28  A  at  case  temperatures  of  100  °C  and  25  °C,  respectively. 

VI.3  Experimental  Results:  UAF-Packaged  Sprayed-Cooled  IGBT 

Table  5  illustrates  the  measured  results  when  the  IGBT  bare  die  is  packaged  according  to  the 
proposed  novel  packaging  methodology.  As  expected,  the  case  temperature  in  Table  3  agrees  well  with 
the  heat  spreader  temperature  since  the  spraying  conditions  are  identical.  The  main  difference  is  in  the 
calculated  junction  temperature  that  is  considerably  smaller  for  the  UAF  power  package  because  it  has 
been  optimized  for  spray  cooling.  The  junction  temperature  for  the  UAF  power  package  is  33  °C,  much 
lower  than  the  59  °C  junction  temperature  of  the  TO-220AB  package  for  identical  power  and  spraying 
conditions.  From  these  results  one  can  envision  that  the  UAF  power  package  should  be  able  to  operate  at 
higher  power  levels  for  the  same  junction  temperature.  Unfortunately,  the  UAF  power  package  was 
render  inoperable  because  of  likely  voltage  spikes  on  the  gate  signals  so  it  could  only  be  tested  to  10.2  A. 


Table  5.  Spray-cooling  data  for  UAF  power  package  at  9.2  kHz. 


■c(A) 

T(heat 

spreader)  (°C) 

T(junction) 

(°C) 

T (fluid  spray) 

(°C) 

PlGBT 

(W) 

0.0 

20.3 

21 

20.0 

0.0 

5.1 

25 

20.1 

13 

10.2 

26.7 

33 

20.8 

27 

28 


Cold  Plate  IR  K3BT 


Spray  Cooled  U  of  A 
Power  Package 

Spray  Cooled  IR  IGBT 


Junction  Temperature  (°C) 


Figure  32.  Collector  current  vs.  junction  temperature. 

Working  with  actual  IGBT  (as  opposed  to  using  thermal  test  heater)  brings  considerable  more  risks 
and  challenges.  The  experimental  set  up  had  significant  amount  of  parasitic  inductances  that  were 
causing  over-voltages  (see  Figure  30).  Measurements  will  be  taken  in  future  prototypes  so  this  problem  is 
not  repeated.  Nevertheless,  this  damage  does  not  take  away  from  the  excellent  performance  of  the  UAF 
power  package.  At  the  writing  of  this  report  another  UAF  package  is  being  developed  for  testing  up  to  its 
maximum  capability  based  on  these  promising  results. 

Figure  32  above  illustrates  in  graphical  form  the  results  of  Tables  3  to  5.  An  interesting  point  is  that 
the  junction  temperature  for  the  UAF  power  package  is  half  of  that  for  the  IRG4PC30KD  package  for  the 
same  current  levels.  This  was  possible  due  to  the  lower  thermal  resistance  value  of  the  UAF  power 
package  built  by  the  researchers  (almost  4  times  lower).  Also,  the  UAF  package  area,  for  the  same  bare 
die,  is  almost  3  times  smaller  than  that  of  the  IR  TO-220AB  package.  These  results  mean  that  potentially 
significant  improvements  in  terms  of  power  density  and  volume  density  could  be  possible  by  employing 
the  novel  packaging  methodology  developed  by  the  UAF  researchers.  The  thermal  advantages  of  spray 
cooling  over  conventional  cold-plate  cooling  are  also  extracted  from  Figure  32. 

For  completeness.  Figure  33  illustrates  typical  measured  gate  signal,  motor  current,  motor  voltage 
and  IGBT  collector-emitter  voltage.  The  electrical  system  was  identical  for  both  cases  with  the  changes 
only  taking  place  in  the  thermal  management  approach 
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Chi  VCE  (100  V/div) 
Ch2  Vge  (20  V/div) 
Ch3  VMotot  (200  V/div) 
Ch4  Ic  (10  A/div) 


Chi  Vce  (100  V/div) 
Ch2  Vge  (10  V/div) 
Ch3  Vm*  (100  V/div) 
Ch4  Ic  (10  A/div) 


Figure  33.  Typical  measured  signals  for  IR IGBT  cooled  by  (a)  Cold  plate  and  (b)  Spray  cooling 
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VII.  CONCLUSIONS 


The  main  objective  was  to  investigate  the  feasibility  of  using  spray  cooling  as  the  thermal 
management  technology  for  power  semiconductor  devices  used  in  power  converter  applications,  in 
particular,  naval  applications.  The  main  conclusions  are  addressed  below  in  relation  to  enabling 
objectives  established  by  the  UAF  researchers  to  accomplish  the  project  main  objective. 

VII. 1  Conclusions  on  Theoretical  Understanding  of  Spray  Cooling  Phenomenon 

Theoretical  understanding  of  the  basics  of  spray  cooling  modeling  involving  phase  change  was  the 
first  enabling  objective  since  previous  research  efforts  have  been  mainly  experimental  with  little  research 
done  on  a  comprehensive  modeling  of  the  phenomenon.  A  comprehensive  modeling  of  the  spray  cooling 
phenomenon  could  be  used  to  design  optimized  spray  cooling  systems. 

As  a  result  of  the  literature  survey  the  UAF  researchers  concluded  that: 

•  The  whole  process  of  modeling  phase  change  under  spray  cooling  is  a  difficult  task. 

•  Modeling  of  the  heat  transfer  mechanism  by  a  single  droplet  on  a  hot  surface  using  a  heat 
transfer  coefficient  and  a  CFD  package  has  limited  application  because  of  the  difficulty  in 
determining  the  heat  transfer  coefficient  under  spray  cooling  from  experimental  work. 

•  Multiphase  flow  modeling  is  a  viable  technique  to  understand  heat  removal  at  the  interfaces 
of  liquid  and  solid  by  a  vapor  bubble . 

•  The  “level  sef  method  is  one  of  the  best  methods  to  solve  multiphase  flow  in  spray-cooling 
systems.  This  conclusion  resulted  from  a  survey  of  the  status  of  computer  modeling  in  spray 
cooling  and  methods  to  solve  multiphase  flow  model.  The  survey  addressed  in  this  report 
was  published  as  a  part  of  the  Selvam  (2005)  paper. 

•  The  modeling  nucleation  and  the  bubble  growth  in  a  thin  film  of  liquid  of  about  70  pm 
should  be  performed  in  a  future  task. 

VII.2  Conclusions  on  Power  Packaging  and  Spray-Cooling  of  Power  Semiconductor  Devices 

The  enabling  objectives  on  this  topic  focused  on  developing  a  power  packaging  methodology 
compatible  with  spray  cooling  since  the  UAF  research  team  felt  that  packaging  must  be  also  considered 
when  applying  spray  cooling  in  order  to  fully  realize  the  benefits  of  this  thermal  management  approach. 
A  technology  demonstrator  consisting  of  an  H-bridge  power  module  driving  a  240- V,  10-hp  dc  motor  was 
assembled  to  demonstrate  the  feasibility  of  the  proposed  ideas. 

Initially,  the  UAF  researchers  used  thermal  test  heaters  on  a  spray-cooling  test  bench  to  evaluate  how 
factors  such  as  pressure,  flow  rate  and  spray  fluid  temperature  affect  the  cooling  of  the  power 
semiconductor  devices.  The  results  demonstrated  that  in  general: 

•  High  nozzle  flow  rates  are  helpful  when  all  other  conditions  are  the  same  because  of  the 
enhanced  conduction  mechanism  for  heat  removal. 

•  High-flow  nozzles  cb  not  always  produce  the  optimum  cooling  efficiency  since  they  may 
hinder  the  phase-change  process  (see  Figure  27). 

•  Droplet  size,  droplet  velocity  and  perhaps  other  factors  clearly  play  a  key  role  in  this  process 
as  it  is  being  analyzed  as  part  of  ONR  Award  No.  N000 14-04- 1-0603. 

The  UAF  researchers  designed  a  spray  box  to  spray  standard  TO-220AB  packages  of  the 
IRG4PC30KD  part  from  International  Rectifier  (having  an  IGBT  and  a  freewheeling  diode)  and  the  H- 
bridge  power  module  packaged  by  the  UAF  packaging  methodology  (see  Figure  23).  The  IRG4PC30KD 
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part  was  selected  simply  because  of  its  availability  in  bare  die  form  (so  it  could  be  packaged  using  the 
methodology  of  Section  III)  and  its  ratings  were  compatible  with  the  electric  motors  available  in  the 
Energy  Conversion  Laboratory  so  experimental  results  could  be  conducted  to  practically  demonstrate  the 
feasibility  of  the  proposed  ideas.  The  IR-packaged  part  was  also  cooled  under  the  same  operating 
conditions  by  a  standard  cold  plate  for  comparison  purposes. 

The  experimental  results  demonstrated  that: 

•  Spray  cooling  is  much  better  than  liquid  cooling  of  a  standard  TO220AB  package. 

•  Simply  spraying  a  conventional  module  does  not  result  in  the  optimum  cooling  efficiency. 
The  package  must  be  also  optimized.  From  an  analysis  of  the  experimental  data,  the 
performance  of  the  spray-cooled  UAF  power  package  far  exceeded  the  performance  of  the 
liquid-  or  even  spray-cooled  conventional  package. 

The  UAF  packaging  methodology  accomplished  the  following: 

•  Wired  bonds,  a  major  source  of  failure  in  power  modules,  were  eliminated  by  replacing 
them  with  solder  balls. 

•  The  package  materials  were  CTE  matched  for  high  reliability. 

•  The  thermal  resistance  between  the  heat  source  and  the  liquid  spray  was  significantly 
reduced  with  the  potential  of  greatly  improving  the  overall  cooling  efficiency. 

•  The  technology  can  be  used  to  implement  chip -scale  packages,  thereby  increasing  the 
system  power  density. 

In  summary,  the  experimental  results  prove  that  spray  cooling  is  an  enabling  thermal  management 
technology  that  should  meet  the  future  heat  load  demands  of  power  semiconductor  devices  used  in  power 
converters.  However,  novel  packaging  techniques  such  as  the  packaging  methodologies  developed  by  the 
UAF  researchers  are  required  to  realize  the  full  potential  of  spray  cooling. 

VIII.  FUTURE  WORK 

Currently,  several  activities  are  underway  as  part  of  a  continuation  project  that  mainly  concentrates 
on  Silicon  Carbide  (SiC)  power  semiconductor  devices  (ONR  Award  No.  N000 14-04- 1-0603). 
Regarding  modeling  the  spray  cooling  phenomenon,  these  activities  are  aimed  at  developing  a  multiphase 
flow  computer  model  using  the  level  set  method.  The  computer  model  will  be  used  to  study  nucleation 
and  boiling  in  a  thin  film  of  fluid  of  about  70-pm  thickness  on  a  hot  surface.  The  UAF  researchers  are 
cooperating  in  the  area  of  spray  cooling  with  researchers  from  the  Air  Force  Research  Laboratory  and 
NASA  Goddard  Space  Flight  Center. 

While  the  results  related  to  cooling  of  power  semiconductor  devices  that  were  achieved  in  this  ONR 
grant  were  a  great  improvement  over  many  other  methods  of  thermal  management,  the  UAF  researches 
feel  that  potential  exists  to  improve  further  the  maximum  heat  flux  that  is  possible  with  this  technology. 
The  current  limitations  are  due  to  the  properties  of  the  used  dielectric  fluid  and  the  fact  that  only  a  small 
portion  of  the  sprayed  fluid  is  converted  to  a  vapor  phase.  Much  of  the  fluid  simply  flows  off  the  surface 
and  contributes  only  to  cooling  by  convection.  On  going  work  is  focusing  on  the  use  of  fluids  with 
greater  heats  of  vaporization  values  as  well  as  methods  to  promote  a  more  complete  vaporization  of  the 
fluid.  An  important  challenge  is  to  be  able  to  practically  realize  the  theoretical  concepts  that  promote 
greater  vaporization  of  the  fluid. 

Current  power  packaging  activities  are  aimed  at  packaging  of  SiC  devices  for  high  voltage  and  high 
temperature.  The  SiC  MOSFET  is  being  targeted  for  high  voltage  packaging  and  spray  cooling  will  be 
used  as  the  thermal  management  approach  since  it  has  a  temperature  limitation  of  approximately  200  °C 
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due  to  the  gate  oxide.  The  SiC  JFET  and  SiC  SIT  are  being  targeted  for  high  temperature  packaging 
(greater  than  250  °C  or  300  °C). 
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